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Abstract 

< 

The quantum dynamics of a two-level system coupled to an Ohmic spin- bath is studied by means of the perturbation 
■ approach based on a unitary transformation. A scattering function is introduced in the transformation to take 
into account quantum fluctuations. By the master equation within the Born approximation, nonequilibrium dynamics 



quantities are calculated. The method works well for the coupling constant < a < a c and a finite bare tunneling 
A. It is found that (i) only at zero temperature with small coupling or moderate one does the spin-spin-bath model 
display identical behavior as the well known spin-boson-bath model; (ii) in comparison with the known results of spin- 
boson-bath model, the coherence- incoherence transition point, which occurs at a c = h[l + r)A/u> c ], is temperature 
independent; (iii) the nonequilibrium correlation function P(t) = (T z (t)), evolves without temperature dependence 
' while (T x (t)) depends on temperature. Both P(t) and {r x (t)) not only satisfy their initial conditions, respectively, 

I ! 

• . and also have correct long time limits. Besides, the Shiba's relation and sum rule are exactly satisfied in the coherent 

. 

' regime for this method. Our results show that increasing temperature does not help the system suppress decoherence 
in the coherent regime, i.e., finite temperature does not favor the coherent dynamics in this regime. Thus, the finite- 
• i-h ' temperature dynamics induced by two kinds of baths spin-bath and boson-bath exhibit distinctly different physics. 

PACS numbers: 72.20.Dp; 05.30.-d; 03.65.Yz. 
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I. INTRODUCTION 

The dynamics of a dissipative two-level system has attracted extensive studies in last decades^, since it can be 
used to describe a large number of different physical and chemical processes. The interaction between the system 
and its environment gives rise to decoherence and dissipation which are also the major stumbling block to quantum 
computation and quantum communication. Generally, there are two kinds of quantum environments, one is boson- 
bath modeled by a set of oscillators (dclocalized modes such as phonons and photons), the other is two-level systems 
(TLSs) reservoir or spin-bath (localized modes such as defects, impurities, nuclear and paramagnetic spins )». Usually, 
people use the spin-boson-bath (SBB) model 

Hsbb = -\&t x +^uib\bi + \^9i{b\+bi)T z , (1) 
i i 

to take into account the system-bath interaction[4-26], where a bosonic heat bath consisting of an infinite number 
harmonic oscillators (denoted by b\' s and b{s) constitutes the environment of a quantum TLS (denoted by the Pauli 
matrices t x and r z ). Experimentally, the spin-bath plays an important role on the decoherence in magnetic cluster and 
semiconductor qubits at low temperature with some interesting features^. For example, GaAs quantum dot electron 
spin qubit, a candidate of solid state quantum computation, loses its quantum memory due to its coupling with the 
surrounding nuclear spin environment being the spin-bath. The abundance of spin-baths in real systems urgently 
necessities an understanding of effects on decoherence. Another simple model for spin-bath is proposed, which is the 
so-called spin-bath composed of an infinite number of TLSs without mutual interaction^— . The Hamiltonian of a 
TLS coupled with a dissipative spin-bath (spin-spin-bath, SSB) reads 

H=-\^T*--Y i ui° l x-\Y,to al * r '- (2) 

I I 

Here t x and r z are Pauli matrices to describe the TLS, a l x and a l z are Pauli matrices for the /-mode of spin-bath^. A 
is the bare tunneling matrix, loi the frequency for the Z-mode of the bath, and gi the coupling constant. The coupling 
between the TLS and its environment is characterized by a spectral density J(w) = ,9; 2( 5( w — w () = 2aui6(uj c — oS) 
with the dimensionlcss coupling strength a, the upper cutoff u> c and the step function 9(x). 

Most researchers are interested in an open system with a small number of degrees of freedom (such as a two-level 
system or an oscillator) in contact with a bath of a complex nature, whose number of degrees of freedom tends to 
infinity. The evolution and its properties of the open system are determined by the coupling to the bath. Do the 
boson-bath and spin-bath have the same effects on the dynamics of TLS? What is the difference between the coherent 
dynamics of SBB and SSB? At first glance, differences of decoherence to the open system could be attributed to 
intrinsic energy level structure of heat baths. The boson-bath may be treated as an infinitely large reservoir of energy, 
since the number of oscillators is infinite and every harmonic oscillator can be excited to the equally distributed states 
without upper limit. On the other hand, the spin-bath may also be treated as an reservoir because the number of spins 
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in the bath is infinite, but there is only a single excited state for every bath spin. Thus, the distinct physics between 
them ascribe to their underlying structures of bath and continue to attract much attention from both theoretical and 
experimental sides^ir— . 

The SSB model was studied by approximate analytical and numerical methods, such as Kubo's cumulant expansion 
method 2 ^, the perturbation theory 2 ^, the resolvent operator approach^, the numerical path integral method 3 -^. The 
main theoretical interest is to understand how the environment, the spin-bath, influences the dynamics of the TLS 
and, in particular, to discuss the common features and the main differences between the dissipative roles played by the 
spin-bath and boson-bath. Some studies show that decoherence is partially suppressed by increasing the temperature 
of spin-bath and temperature plays, though weakly, a positive role in maintaining coherent dynamic a 30 ' 31 , which 
is in contrast with the conclusion of the cumulant expansion method that the decay rate of TLS is temperature 
independent^. 

For instance, the study of noninteracting blip approximation (NIBA) can give the population difference 
P^nib a 31 ^ 8 ■ Its Laplace transform P(s) for spin-bath is given 

P( S ) = [ S + f(s)}-\ (3) 

where 

poo 

f(s) = A 2 / cos(Qi(t)) exp(-st - Q 2 (t))dt, (4) 
Jo 

and the function Q\ and Qi are given by the relations 

Qi = sin(arf)tanh(— J-^dw, (5) 

Q 2 = /~[l-cosM)]^dw. (6) 
Jo 

As shown in Ref— , their simulations as well as the solution of the NIBA equations, indicated that the diffuse coherent- 
incoherent boundary shifts to stronger coupling as the temperature is raised, which is contrary to the known behaviors 
for boson-bath(With increasing temperature, a c decreases quickly). It stands a striking contrast against the general 
belief that the heat bath with infinite degrees of freedom (thermal reservoir) leads to the dissipation and decoherence 
of the open system and the increase of temperature quenches or docs not favor the coherence. Despite of some different 
arguments for the coherent-incoherent transition in this model, these works highlight that the decoherence of open 
system depends crucially on the underlying nature of spin bath. Thus we would reconsider the topic by an analytical 
method based on a unitary transformation. 

In this work we present an analytical approach for calculating the dissipative quantum dynamics of the SSB model. 
It works well for the coupling constant < a < a c and a finite bare tunneling A, and could reproduce nearly all 
exact results obtained by various analytical and numerical method a 23 ' 39 . It could explain the crossover between the 
coherent oscillation and the incoherent behaviors and also allow us to resolve the controversial claims in the literature. 
Throughout this paper we set H = 1 and ks = 1. 
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II. UNITARY TRANSFORMATION 

A unitary transformation, which is defined as H' — cxp(S)H exp(— S), is applied to H and its aim is to take 
into account the correlation between the TLS and its bath. To this end, the following form for the generator is 
propose d 23 ' 39 , 

Here a /-dependent function is introduced in S and its form will be determined later. 
The transformation can be done to the end and the result is 

H' = cxp(5)Hcxp(-S*) 

= -^Ar g cosh^ vtia l v ] - - Air y sinh[]P viia l y ] 
i i 

~2 ^2^1 cos M + 9i sin (wz)]cr' + ir 2 Y^i^i sin ( u ') - 9i cos(vi)]a l x , (8) 
l i 

where vi = gi^i/ui. The transformed Hamiltonian H' may be divided into three parts, 

H = H' + H[ + H' 2 , 

H' = - l - V Ar x - l -Y.^ ( 9 ) 

i 

with a renormalized frequency for the bath, 

lu[ = u>i cos(w/) + gi sin(tv), (10) 

and a renormalized factor for tunneling 

r) = Tt b (pb cosh[^2 viicr l y ]) = cxp(^ ln[cos(i>/)]). (11) 



Here pb = exp(— j3 H g) /Ti exp(— (3Hb) is the equilibrium density operator of bath spins (i/g = ~\ X)/ ^j !) an d Tr^ 
is the trace operation with respect to the bath. is the unperturbed part of H' and, obviously, it can be solved 
exactly. The eigenstate of H' is a direct product: |s)|{±l;}), where \s) is the cigenstate of t x : \si) = -j- ( J ^ or 

\s2) = ( \ ^' anc ^ I "l =1= 1 ^ } ) i s the eigenstate of bath spins: {±1/} means that the eigenvalue of a\ is +1 or — 1. 
In particular, |{+1}) is the vacuum state with eigenvalue of every a z being +1. Then, the ground state of H' is 

\go) = \si)\{+i}). 

The other terms in H' are 

H 'i = \ T * s ' m M - 91 cos(u z )] - ^AiTy sm{vi)icr y , (12) 

l i 

H' 2 = -\&t x fcosh[^ «jioi] - r\ 



(13) 



H[ and H' 2 are treated as perturbation and they should be as small as possible. For this purpose 77 is determined in 
Eq.(fTTj) to make TxsiPBH^) = 0. Besides, £2 is determined as 

Vl = iL^ = tan -i — j—r, (14) 
uji lui + r)A 

and because of this definition one obtains 



Here we use the following definition for spin operators: a l + = \{cr l x — i<J l y ), a l _ = \{o l x + i<r l y ), t + = \{r z + ir y ), 
t_ = \ij z — iT y ), since the ground state of Zth bath spin is a\ = +1 but its excited state is a\ = — 1. It is easy to 
check that H[\go) = because H[ contains the rotating- wave terms only. The counter- rotating- wave terms in original 
Hamiltonian contribute to the renormalization in H' and its form in Eq. (12) disappears because of the functional 
form of £1 in Eq. (fT4)) . At the same time, the bare coupling gi in H has changed to an effective coupling in H[. These 
are essential in our approach. 

Substituting (fT4)) into (|T0| and (fTTj) . one obtains the renormalized frequency 

, _ uji{uji +i]A) +gf 



and the renormalized factor of tunneling 



r] = exp 



-IE- 
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(16) 



(17) 



In our treatment H' is the unperturbed Hamiltonian with the renormalized parameters. H[ is the perturbation which 
contains flip (transition) of single bath spin. H' 2 will be neglected because it involves multi-flip (multi-transition) of two 
or more bath spins and its contribution to physical quantities is 0(gf) and higher. Thus, in the following treatment, 
H' Ki Hq + H[. 

The renormalized factor of tunneling 77 is important to determine the physical property of the coupling system. In 
order to make the summation over I of physical quantities, without loss of generality, a constant density of state is 
introduced, p(u>i) = po: 

J2f("i,9i)= / du lPo f(ui,gf), (18) 
1 J o 

gf = 2aui/po, (19) 

where f(uii,gf) is any function of lui and gf. It can be checked that the above treatment is in agreement with the 
Ohmic spectral density J(w). In following calculations the dimcnsionless quantity pouj c is treated as the total number 
of bath spins which goes to infinity in the thermodynamic limit. For example, the summation in Eq. ([17[) can proceed 



as 



dx \ _ / cj c + r?A 

(x + rjA/uc) 2 ) ~° XP V ??A ' u; c + riA 



77 = exp ( -a / — ; - — — ) = exp [ -a In — — h - — ) , (20) 



where x = loi/<jj c . In the scaling limit A -C lo c one can get an explicit solution 77 = (eA/uc)"^^ 1 *^ which leads to a 
localization point at a = 1, 77 = for a > 1. As the coupling increases, rj decreases smoothly to zero. 
Eq. (j2"0)) is exactly the same as that of the SBB model (Ohmic bath) for ground state at T = 0^2, 



r] B = exp 



-V g ' — , coth(^)) . (21) 
2^(cj / +? ? A) 2 y 2T' I 



rj of the SSB model is temperature independent but rj B is temperature dependent. That is because every bath spin 
has only a single excited state but every oscillator in the boson-bath has infinite excited states. In addition, it is seen 
that t\b decreases with the increase of temperature. In other words, the raise of temperature of bosonic bath leads to 
the faster loss of coherence. 

As the transformation in Eq.® has been done without approximation, one can calculate the upper bound of the 
ground state energy of the coupling system by 



E g /u c = -\]A/u c - - Y^'l - ui]/u c 

1 

1 . . a f 1 , x 2 + 2xrjA/uj c 1 A , a 

= -^Aju] c - 2 y o dx 



(x + V A/u c ) 2 ' 2 VA/UJC 2 1 + ?? AK' (22) 

Note that the ground state of H is exp(— s)\go), it is the ground state of interacting system. The effect of fluctuating 
environment has been taken into account in the treatment. 

III. NON-EQUILIBRIUM CORRELATION FUNCTION 

The density operator in Schrddinger representation is psB(t) with Hamiltonian H, where the subscript SB indicates 
that it is the density operator for the coupled two-level system and bath. For transformed Hamiltonian H' the density 
operator is p'g B {t) = e s psB{t)e~ s . We treat H' as the unperturbed Hamiltonian (Eq.(5)) and the density operator 
in the interaction picture is 

p'£ B (t) = exp(iHl,t)p' SB {t) expi-iH^t). (23) 

The equation of motion for p'g B (t) is^ 

f t p'£ B (t) = -i[H' 1 (t),p'i B (t)]. (24) 
H[{t) is the perturbation H[ (Eq.(15)) in the interaction picture, 

H[(t) = -J2 V l [e-' l( ^- nA)t T + <7 l _ + e^-^T-a 1 ^ . (25) 



where Vi = r)Agi/[(u>i + 77A) 2 + gf] 1 ^ 2 - Our procedure for solving the equation is to write^i 

PSB(t) = Ps(t)pB+p' c I (t), (26) 

where p'g(t) = Tr B p'g B (t) is the reduced density operator for TLS. Here p'J(t) is the cross term of system and bath 
operators representing the correlation because of their interaction, which is the order of gi and higher. TibPc (t) = 0. 
Now Eq. (|!M|) can be integrated as 

ti®PB-ti{G)PB + £{t)-pf!{G) 

= -i f ''dt'[H[(t'),p'i(t')p B ]-i f dt'[H[{t'), P '!{t')]. 
Jo Jo 

After applying the Ttb operation to both sides, we have 



P'i(t) -pf£(0) = -* f dt'^ B [H[{t'),p'^t')]. (28) 
Jo 

Other terms in Eq.(27) are 

p' c I (t)=p' c I (0)-i f dt'lH'S'lp'^pB]. (29) 
Jo 

Substituting this equation into Eq.(28) we get the master equation for p's(i)r^ 

d p'i(t) = -iTrslHUtlp'^O)) f^BlH^tllH^p'^pBldt'. (30) 



(27) 



o 



dt 

where all higher order (than gf) terms are neglected. 

At t = 0, the usual initial density operator is psb{0) — ^ q q J Pb- Then one can get the initial condition for the 
calculations: p'g B (0) = Pg B (0) = e s psB(0)e~ s leads to 



p'i(0) = Tr B p' SB (0) = ( I o ) , tf(0) =[s,(H)pb 



(31) 



where we stop at the first order of gi (Eq.([7]) for S). 

The details of perturbation calculation are listed in Appendix. Note that temperature dependence (tanh(^)) ap- 
pears explicitly in the derivation. We obtain that the diagonal elements in the reduced density matrix are independent 
of temperature, nondiagonal elements are dependent on temperature. The solution of the reduced density operator 

*<"-(£&)'■ 

1 f' 00 1 



M) - fat) = ^ J ^du> 



l u)-u)[+iO+ 

1 



cu — 7]A — V, 

(32) 



The real and imaginary parts of J^i Vf/fa — *0 + — w arc denoted as R(ui) and 7(0;), respectively. They are 



^/{uji + 77A) 2 + gf uj^/iuji + 77A) 2 + gf - wj (wj + ?/A) - g\ 
(77A) 2 / cj c w ^ [|w|(o; c + 77A)" t 



= -2a w . ^ ^— hi 

+ 7/A 0J C + 7/A 07 + 77 A 



77A(w c — Oj) 



7Tr, 2 A 2 gf w,(w,+J7A)+flf 



(?7A) 2 , , 

= 2cnra; . w 34 
(w + 77A) 2 v y 

The decay rate 7(0;) is not dependent on temperature but on frequency which agrees with the conclusion of 
Nitzan and Silbey's paper—. Its nonmonotonic behavior exhibits a maxima at uj = 77A which is distinguished from 
the rotating- wave approximation (RWA) results (jrwa = 2anuj). The non-equilibrium correlation function P(t) is 
defined as 

P(t) = (T z (t)) = Ty s Ti- B (psB(t)T z ) = Tr s Tr B (p' SB (t)a z ) 

= Tr s Tr B {[p' s (t)p B + e^'p^e"^) = Tr s (p' s {t)r z ) 

1 r \ t , s 1 f u " , 7(w) cos(wt) , „ 

= Pii(t) - p' 22 {t) = - / duo- n ' ' 2 , v 35 

since Tr^ps = 1 and Ttbp' c (f) = 0. Here {r z (t)) is used to denote the average TrsTrs(psB(i)T z ), which is the 
population difference. The last equality in Eq. (|35| comes from Eqs. (f32jl with the use of the Kramers-Kronig relation. 

Note that P(t) is temperature-independent, which means that the increase of temperature does not favor nor 
suppress the coherent dynamics, which is distinguished from the temperature effects in the SBB model(the coherence 
loses quickly with increasing the temperature). Fig. 1 shows the P(t) versus At relations with A/oj c = 0.1 for different 
couplings a. From the evolution behavior of the SSB, it is found that, with increasing coupling, the dynamics exhibits 
from the damped coherent-oscillation for the weak coupling to incoherent decay for moderate coupling. Fig. 2 shows 
the P(t) versus r]At relations with a = 0.1 and different tunneling. The curves show a similar scaling behavior for 
different tunneling A. 

The integration in Eq. (|35| can be done approximately by the residue theorem, 



P{t) = cos(w i) exp(-Tt), (36) 

where ujq is the solution of equation 

uj - 7/A - R(u ) = 0, (37) 
and 7 is the Wigner-Weisskopf approximation (WWA) of 7(0;): 

7 = 7(07 = 77A) = — anr/A. (38) 
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Notice that ujq — A = (77 — 1)A + R(u>o) is the level shift induced by the spin-bath coupling. In weak coupling case, 
the integral function in Eq. (|32p of P(t) possesses two complex poles which result in damped oscillation dynamics^. 
The real part ujq represents the frequency of coherent tunneling. With increasing coupling, ujq become smaller. The 
solution ujq of Eq. ([57|) is real (u/q > 0) only when a < a c , and 



is determined by ujq — in Eg. (1371) . It becomes the well-known result a c = 1/2 for the SSB in the scaling limit 
A/aj c <C 1. For a > a c there is no real solution ujq and it means that a = a c determines the critical point for a 
coherent-incoherent transition in contrast to the diffuse boundary region between the coherent phase and incoherent 
one predicted by numerical path integral treatment^!. Besides, the critical coupling is independent of temperature, 
so there is no boundary shift as the temperature is raised, which is in contrast with the known behavior of an Ohmic 
bath of bosons and the NIBA results with weak-temperature dependences^. Thus, it is found that the increase of 
spin-bath temperature does not favor the coherence. 

In Fig. 3, a phase diagram shows the relation between temperature and coupling. In comparison with numerical path 
integral results, and the NIBA results of both the SSB and SBB models, our analytical results show the crossover from 
the coherent oscillations to incoherent decays is independent of temperature. The dashed line shows the coherent- 
incoherent boundary predicted by the NIBA with an effective harmonic bath of temperature-dependent spectral 
density for ui c = 20 A (All data are taken from Ref£L). Thus, the NIBA solution for the spin-bath indicates that the 
coherent-incoherent boundary shifts to stronger coupling as the temperature is raised, which might be questionable 
since NIBA is not reliable for small values of cutoff frequency and low temperature regim o 1 ' 2 ' 31 . From Fig. 3, one 
can see that at lower temperature, the critical coupling for the SSB obtained by the NIBA is about 0.75 which is 
larger than a c — 0.5 for the SBB model for A/uj c <C 1 at T = 0. However, at zero temperature, from the form 
of bath correlation functions of the two models by the NIBA shown in the followin g 31 i 38 , the SSB and SBB model 
yield identical results which is the same as the conclusion obtained by the resolvent operator approach that there is 
the same underlying physics for both two models at zero temperature 3 ^. Additionally, the simulations by numerical 
path integral is also shown for comparison in the shaded area. The parameter space is obtained by an experiential 
method that P(t) has a small negative lobe that does not fall below —0.01 which is stated in Ref. f 3 -^]. Consequently, at 
moderate or higher temperature, the width of this area is considerably broad whose parameters is corresponding to the 
criterion 3 -^. In contrast, our approach formulates P(t) for the SSB model obviously without temperature dependence. 
In the scaling limit, one can readily get a c = 0.5 which is consistent with the exact results^i 3 ^. As a consequence, our 
result turns out that the coherence of the TLS does not benefit from the increase of temperature. 

Since Eq. (|35p is temperature independent, these conclusions hold also true for the finite temperature. This is totally 
different from those of SBB, as the coherent oscillation of SBB disappears quickly with increasing temperature. On the 
other hand, NIBA results for boson-bath indicate that the transition temperature drops very quickly as the dissipation 




(39) 
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increases which is shown in Fig. 3. The P(t) of the SBB model is also calculated by our approac h 23 ! 39 , 

t Jo i w - V A - Rb{u)\ z + T B \ U ) 
where and 7b (w) are the real and imaginary part of VJ 2 coth(^)/(w — ojk — i0 + ), respectively (For boson 

bath, Vi = t]b Agi/ (() B A + w/)). They arc given by 

r b (u>) = -2* {iibA)2 J " c a - - ln r n%+^) 1\ 

+ V-^ , (41) 

lo - ui exp(w fc /T) - 1 

7b(w) = 2q^ ( ^ A) " coth(^), for < u < u> c . (42) 

In order to show our method explicit, we compare the result of our approach to the results on the spin-boson model at 
T = from numerical renormalization group theory^, which is shown in Fig. 4. It is found that from weak coupling 
to moderate coupling, our result is in good agreement with those of the numerical method. However, there appears 
a difference near the coherent-incoherent transition. 

At T = 0, the population difference of the SSB model is the same as that of SBB model. It confirms that they 
exhibit the same dynamics at zero temperature which coincide with the findings obtained by Shao and Hanggi 3 ^. 
The decay rate j B (u)) is temperature dependent which is consistent with the known results^^. By comparing 7(0;) 
in Eg . ([35)) (without temperature factor) at finite temperature with 75 (w) in Eq. (|40"]) fwith coth(^y)), we find that the 
difference of dissipative roles between the two kinds of baths ascribes to the two distinct level structures for each bath 
degree of freedom and available states distribution. 

The difference between the two models can be traced back to the restriction of the thermal induced excitation 
possibilities of any bath degrees of freedom. For spin-bath, there is only a single level of thermal excitation in each 
individual two-level system of bath, while for boson-bath, there are infinite levels in each individual oscillating mode. 
In the boson-bath model, thermal excitation of many levels of a single bath degree of freedom is one of important 
mechanism for decoherence, whose excitation number can be represented by + 1 = coth(wfc/2T), while the spin- 
bath model lacks in this structure. Thus, our approach formulates the decoherence measure P(t) for the SSB model 
obviously without temperature dependence. 

( T x(t)} can be calculated in a similar way as Eq. (|35l) . 



< r x {t) >= TT S Ti-B{p' SB {t)e S T x e- s ) 

= TrsTrs ([p' S (t)p B + e' 1 ^ p 1 ^ {t)e lH '^][r x cosh(X) + ir y sinh(X)]) 
= Tr 5 Tr B (p' s {t)p B T x cosh(X) + ip'/^e^V^e"^* sinh(X(t))) , 

where X = J2i v i( (jl - ~ ° + ) an d X(t) = Yli Vi{a_e~ u ° lt — <j l + e' lullt ). The trace operation related to the cross term 
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p'J (i) can be done with Eq.(Al6) and the result is 

(r x (t)) = r?tanh(|^) {1 - cxp(-2 7 t)} - I ^viVi tanh(^) sin( W/ <) sin^Ai) 

z 

+1 Y ^ 2 [tanh(^) - tanh(^)] 1 = ^ " ^ )t] 
A ^ 1 1 y 2T' y 2T n uji - rjA 

, 1 t/2 t u^ A \ ~ ^A)[exp(-27i) - cos[(w ; - r/A)i]] + 27sin[(£j; - ??A)i] 
+ tanh( — ) (43) 

where we have taken into account the terms up to the second order of gi. Here, the temperature plays some role, and 
Eq. ((33]) leads to correct long time limit rytanh(^). One can easily check that the initial conditions P(0) = 1, and 
( T x(t = 0)) =0 arc well satisfied. Besides, 

r?A 

(r x (t -> oo)) = 7? tanh(i^), P(t -> oo) = 0, (44) 

which are the correct results for thermodynamic equilibrium state. 

In our work, the dynamic behavior of the central spin is mainly determined by the real and imaginary parts of the 
self energy J^i I ( w — *0 + — cj ; '), R(iS) and 7(0;). Eqs.(33) and (34) show that both R(u>) and 7(0;) are temperature 
independent and this leads to our conclusion that finite temperature does not favor the coherence. However, the 
NIBA may lead to the population difference Pn 1 BA^t) (Refs^i and^S) with the Laplace transform P(s) 

P(s) = [s + f(s)}- 1 , 

where 

f{s) = A 2 / cos(Qi(t)) exp(-si - Q 2 (t))dt, 



f ■ , N , , w , JUS) , 

Q\ = / sm(wt)tanh(— ) — 
Jo 2T w 







1 - COs(cjt) ]—hr-d(jJ. 



Nevertheless for bosonic bath, the two functions become 



sm ( uit ) — ^— did , 



Q 2 = / [l-cos(cjt)]coth(^)^da;. 
Jo 2r w 



Then, Ref.28 claimed that for the spin bath one may introduce an effective spectrum J e ff(uj) = J(w) tanh(cj/(2T)) 
and the coherent-incoherent boundary shifts to stronger coupling (larger a) as the temperature is raised, which is 
contrary to the effect of bosonic bath (with increasing temperature, a c decreases quickly). 

In Refi^S, a polaronic transform is used and, then, the resolvent operator approach is applied to the transformed 
Hamiltonian. It is founded that for the SSB model P(t) is effectively temperature independent at low and high 
temperature, which is the same as ours. For finite temperature, with the similar formulation as those of NIBA Ref^ 
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concludes that temperature plays, though weakly, a positive role in maintaining coherent dynamics, which is different 
from our result of temperature independence. 

Generally speaking, the polaronic transformation (Rcfs^S and^) and the NIBA (Refi^i) lead to the second order 
perturbation with perturbation parameter A, since /(s) is explicitly proportional to A 2 and Q\ and Qi are A 
independent. However, based on the transformed Hamiltonian our approach is the second order perturbation with 
renormalized coupling parameter Vi which has included the effects of the renormalized tunneling arising from the 
coupling to the bath. For weak coupling, Vi — gi- In calculations we take into account all the second order terms 
of gf but neglect higher order terms included in H' 2l which are the multi- flip-flop processes of bath spins of different 
modes, that is, two or more bath spins flip of different modes at the same time. Thus, we provide some general checks 
based on some known results. The treatment can be justified by the shiba relation and sum rule. 

IV. CORRELATION FUNCTION AND SHIBA RELATION 

Since exp(S')r 2 exp(— S) = t z , the retarded Green's function can be written as 

G(t) = -i6(t)Z- lr k {cxp(-PH'){cxp(iH't)T z exp(-iH't),r z }} , (45) 

where H' w H' + H[ is the transformed Hamiltonian. The usual notation for the Fourier transform of G(t) is 
G(u) = ((t z ;t z )), where ((A;B)) denotes the retarded Green's function for operators A and B which satisfies the 
following equation of motion, 

uj((A;B)) = ({A,B}y+(([A,H'};B}), 
({A, B}}' = Z-'Tv {exp(-/?tf'){A B}} . 



So the solution for G{uS) is 



GH u-r,A- £ z V?/{w - w {) + u + vA-Ei V?/{u + "0 ' (46) 



Thus, the symmetrized correlation function 



C(t) = ^Tr {exp(-j3H)[T z (t)T z +T z T z (t)]}/Z 



1 

'2^ 



rfwImG(w) exp(— iut) 



1 f°° 7 (w) 

= — / duj- — — — cosfwt). (47) 

7T J [w-77A-i?(w)] 2 + 7 2 (w) 

The susceptibility at zero temperature is related to the Green's function by x"(oj) = — ImG(u;)Sgn(w) as follows, 

X 1 ; [w- A ? /-i?(w)] 2 + 7 2 H [u + A ?/ + i?(-w)] 2 + 7 2 (-w)' V ' 

The static susceptibility \q can be extracted with a Kramer's-Kronig relation and a function-dissipation theorem 

Xo = ~f<*XH (49) 
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One can check the Shiba-relation^i^ 



CM 




(50) 



w->o J(w) 



Note that our normalization condition duiC(ui) = 1, i.e. the sum rule is another check for our approach. We 



provide results for the Shiba-relation and sum rule in the Ohmic case for various values a and A in Table 1. It 
turns out that the Shiba-relation is exact satisfied in numerical precision in the coherent regime. Outside the regime, 
the agreement is still good but no longer exact. Approximations schemes like NIB A or numerical methods based on 
Monte Carlo cannot be used to verify the Shiba relation since they fail to predict the correct long-time behavior—. 



The physics of the SSB model is studied by means of the perturbation approach based on a unitary transformation. 
Analytical results of the quantum dynamics, described by the reduced density operator p(t), is obtained for both 
the scaling limit A/w c <C 1 and the general finite A/uj c case. P(t) is temperature independent while < r x (t) > 
is temperature dependent. Moreover, the decay rate is temperature independent which is in good agreement with 
the conclusion of Nitzan and Silbcy's paper—. It is found that the transition from coherent to incoherent dynamics 
happens at a c = ^[1 + t]A/ui c ], which is temperature independent. Our results have answered the two problems 
mentioned in the introduction. Even though the SSB model has the same dynamics as the SBB model at T = in 
the coherent regime, they displays distinctive dynamics at T > 0. Furthermore, in the boson-bath, the population 
difference decreases fast with increasing temperature, while it is independent of temperature in the spin-bath. Besides, 
the dynamical properties obtained by our approach can both well satisfy initial conditions and reasonably obtain the 
thermo dynamical limits. The conclusion that the coherent oscillation (the population difference) does not depend 
the temperature for the SSB model is not a bad news to the study of quantum information processing in the low 
temperature regime by nanomagnets and nuclear spins. 

If P(t) is calculated by the polaron transformation, the second perturbation theory in the tunneling matrix element 
is applied. Evaluations about P(t) reproduces the same expression of the NIB A which is seen in Ref.£2£](H. Dekker, 
Phys. Rev. A35, 1436(1987)). Note that the NIBA self-energy function is the second order A because there only 
exists a prefactor A 2 in Eq.(4), both Q\ and Q2 are independent of A. However, our approach is the combined second 
order perturbation with renormalizcd coupling Vi which has included the effects of tunneling A and coupling to the 
bath a. On one hand, the contribution of higher order A has been taken into account due to the form of with 
renormalized tunneling. Physically, the self energy function Vj 2 /{ui — i0 + — w ; ') with the renormalizicd coupling 
Vi has been included the higher order of A 2 , whose real part and imaginary one exhibit it explicitly. On the other 
hand, the coupling to bath leads to diagonal transitions and off-diagonal transitions. All diagonal transitions have 
been accounted by 77 in Hq, while non-diagonal transitions have been considered to the first order Vi which has been 
included into H[ and give effects of order a. Obviously, if we make = 1 for any k , the same conclusion as the 



V. 



SUMMARY AND DISCUSSION 
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NIBA can be drawn. Since our transformation is different from the usual polaron transformation, we come to different 
conclusions. 

In our treatment two approximations arc applied. One is the omission of the perturbation term H 2 corresponding to 
multi-phonon non-diagonal transitions such as b k b k , in the spin-boson bath model or multi-bath-spin flips on different 
modes at the same time such as ' a^, in the spin-spin-bath model due to its contribution to physical quantities 
O(g^) or higher orders. This approximation can be justified by the Shiba relation and sum rule (or the initial 
condition P(t = 0) = 1). The other approximation is usual Born approximation for deriving the master equation. 
The generalized Shiba relation and sum rule provide a good test of our method, not just at low frequency, but at all 
energy scales and for values of a beyond weak coupling and a finite A. 

The quantum dynamics of both the SSB model and SBB model are studied analytically by the perturbation method 
based on a unitary transformation. Our approach can be justified by our numerical results: (l)The long-time limits 
of P(t) and {r x {t)) are correct, namely, the expectation value (r z (oo)) vanishes as expected and (t x (oo)) goes to 
the value of thermodynamical equilibrium. Besides, the initial conditions are correct. (2)The coherent-incoherent 
transition point a c = 0.5 is the same as the known result of previous authors. (3)The shiba's relations of both SBB 
and SSB models with Ohmic spectrum have been checked well within a wide parameter range (to see Table I and our 
previous work Ref^). 

Here arc a few words about the key ingredient of the approach. The key point of our treatment is the unitary 
transformation with generator Eq.(7), where a parameter is introduced. After the transformation a perturbation 
expansion has been performed. If = for any k, that is, without the transformation, the perturbation expansion 
would be similar to the standard weak-coupling expansion (Bloch-Rcdficld theory). In addition, if = 1 for any 
k, then our transformation is the usual polaronic transformation and the perturbation expansion is for the small 
parameter A which is equivalent to the NIBA (H.Dekker, PRA35, 1436(1987)). Our choice for < £ k < 1 (Eq.(14)) 
is between them and thus is an improvement on the analytical methods. 

The purpose of our unitary transformation is to find a better way to divide the transformed Hamiltonian into 
unperturbed part H' , which can be treated exactly, and perturbation ones H[ + H' 2 , which may be treated by pertur- 
bation theory. If one treats the coupling term in the original Hamiltonian H as the perturbation, the dimensionless 
expanding parameter is gf/u/f. For Ohmic bath s = 1 it is 2a/ u> which is logarithmic divergent in the infrared limit. 
By choosing the form of r\ (Ea.(|20p) and introducing the function in the unitary transformation it is possible to 
treat H[ and H 2 as perturbation because of the following reason. On account of the form of rj in Eq. ([2U)) H' 2 can 
be treated as perturbation because its contribution is zero at second order of gi. The effect of the coupling term 
in H' (H[) can be safely treated by perturbation theory because the infrared divergence in the original perturba- 
tion treatment for H is eliminated by making choice of the function form The expanding parameter (s = 1) is 
gf^f/ujf ~ 2auj / (uj + r/A) 2 , which is finite in the infrared limit. This approach works well for the low-temperature 
coherent region and the tunneling < A < lu c . It is quite tractable and physically clear, it produces nearly all results 
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which agree with exact ones obtained by various complicated methods in the SSB and SBB model. Thus it may be 
easily extended to more complicated coupling systems. 
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In this Appendix we list the details of solving the master equation (|30p. The first term at right side of Eq. (|3"0"|) is 



and 10547126. 



Appendix 





(Al) 



The integration in Eq. (j3"0"|) can be done as follows, 




- [1 + tanh(^)] [r-p'i(t')r + - p'^t')r + T.]\ exp[z( W; ' - „A)(* t')] 
+ [[l+tanh(M)][r + r_p'/(t')-r_p'/(t')r + ] 

- [1 - tanh(M)][r + p'i(t')r- - p'£ (t')r-r + ]] e Xp [-i(^ - V A)(t t') 



(A2) 



Thus, Eq.([30|) can be solved by the Laplace transformation, 





0) - T+p' s (p)r. 



_ [i + tanh(^)][r_p'i(p)r + - //j( P )r + T_] 
+ [[1 + tanh(^)][r + r_p'i( P ) - r_///(p)r + ] 



p - - r)A) 



1 




(A3) 
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If we denote 



the solution of Eq. (A3) is 



i 

* 

p 



p'A-p& = 



1/2 



1/2 



p— — rjA) 



1/2 



1/2 



p + E 



i p-M(w! — t;A) 



p + E 



E^ 2 t a nh(M)^___ 



1 



^AF 



(A4) 
(A5) 

(A6) 
(A7) 



Using the relation between Schroedinger and interaction representation (23) and making the Laplace inverse- 
transformation, we can get 



Pn(t)+P22(t)=Pll(t)+P22(t) = h 

Pn(t) - PM = coa(riAt)(j>^(t) - pg(t)) - 1 sinfoAi) (t) - p' 2 \(t)) 



(A8) 



p + ir/A + J2 



I p+iuj> 



p - ir]A + 



I p—iuj^ 



(A9) 



Pi 2 [t) ~ p' 2 i(t) = cos( ?7 At)(p' 1 / 2 (t) - p' 2 {(t)) - isiii(TjAt)(p&(t) - p&(t)) 

1 1 



1 



p + irjA + 



p — irjA + 



v? 

I p—iuj^ 



P'l2(t) + P2l(t) = Pl2(t) + P2l(t) 

e pt dp 



1 

2~7Ti 



(A10) 



(All) 



V (l + 2 Ez p i+(j;%A)i) 

The integration path is on a line parallel to the imaginary axis of complex p plane from p = + — too to p = + + ioo. 
The integration in (All) can be re- written as 

[ 1 , 1 

p+i(uj — rjA) p— J?A) 



2^ / 



E ; ^tanh(^)^ 



1 

2ttI 



3 pt dp- 



E; ^ tanh( M 


[ 1 1 1 


} 

1 


p+i(u)[— 7jA) p— ?jA) 
>|W,'+t;A [ 1 


' 2r^A p-H(wi — ??A) p— i(u/ — ^A) 




E^ 2 


1 1 1 




p+i(w ; '-jjA) p—i(wj —jyA) 



P P + Ei ^Z 2 p+ l (J-r,A) + p-iK- 



(A12) 
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The pole point of the first term in {...} is with the residue tanh(?/A/2T). The pole point of the second term is a real 
number, which can be determined by letting l/\p+i{u>i— r)A)] + l/\p— i(u>',— r?A)] = 2p/[p 2 + (u>' l — rjA) 2 ] = 2tt8{uj' 1 — r\ A) 
in summation. Thus, 



r?A 

+ P21M = tanh(^) (1 - exp(-2 7 t)) ; 



(A13) 



where 



— airrjA. 



(A14) 



The integration in (A8) and (A9) is more complicated than that of (A10), because the pole points in complex p 
plane have non-zero real and imaginary part. We change the integration variable in (A9) and (A10) from p to lu: 
p = + + iu) = i(uj — i0 + ), and the result of (A9) is listed in Eq. (f52")l . In order to perform the calculation in Eq.(29) 
the expression for p'g(t) is should be given by 



P5W = ^{l + (l-^ 27t )tanh(^)r : 
7(w)c?u; 



[lj- V A-R(oj)} 2 +j 2 (uj) 



g-iC^-^A)* 



e i(u-riA)t T _ 



Then, the cross term is 

p'J(t) = 



' I 



(A15) 



1 



(A16) 
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Figure Captions 

Fig. 1 P(t) as a function of At for A/uj c = 0.1 and different tunneling a = 0.05 (solid line), 0.1 (dashed line), 
and 0.25 (dotted line). 

Fig. 2 P(t) as a function of rjAt for a = 0.1 and different tunneling A/ui c = 0.01 (solid line), 0.05 (dashed line), 
0.1 (dotted line) and 0.2 (dashed-dotted line). 

Fig. 3 Phase diagram for showing the relation between temperature and coupling for the SSB model in comparison 
with the results of SSB and SBB models by the NIBA. The shaded area is those results by numerical path integral 
method, which indicates that the coherent-incoherent boundary is diffuse. The solid line is our result obtained by 
our method. The NIBA results for spin bath is shown in the dashed line and those for boson-bath in dashed dotted 
line. In addition, the inset displays the coherent-incoherent boundary with relation between a c and A/io c for the SSB 
model. 

Fig. 4 S z (t) = P(t)/2 as a function of ui c t for A/tu c = 0.1 and different coupling a = 0.05 (solid line), 0.1 (dashed 
line), 0.25 (dotted line) and a c = 0.5121 (double-dotted dashed line). The data of numerical renormalization group 
are also shown for comparison in red smooth lines. 
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Table Captions 



TABLE I: Representative results from the numerical solution with parameters chosen by the spectral density 
J(lo) = 2aijj<d(uj c — ui) and with the controlling precision 10 -5 for iteration. R = [lim w _i.o C (ui) / J (uj)] / (xo / 2) 2 ■ The 
numeric error for the Shiba relation and sum rule is at least less than 10 -6 and can be improved by increasing the 
accuracy of numerical calculations. 
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